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Abstract 

We present the first results of a self-consistent solution of the semiclassical Einstein 
field equations corresponding to a Lorentzian wormhole coupled to a quantum scalar field. 
The specific solution presented here represents a wormhole connecting two asymptotically 
spatially flat regions. In general, the diameter of the wormhole throat, in units of the 
Planck length, can be arbitrarily large, depending on the values of the scalar coupling 
^ and the boundary values for the shape and redshift functions. In all cases we have 
considered, there is a fine-structure in the form of Planck-scale oscillations or ripples 
superimposed on the solutions. 
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Wormholes are topological handles in spacetime linking widely separated regions of a single 
nniverse or “bridges” joining two different spacetimes. Interest in these conhgnrations dates 
back at least as far as 1916 |]^ with pnnctnated revivals of activity following both the classic 
work of Einstein and Rosen in 1935 0 and the later series of works initiated by Wheeler 
in 1955 p|. More recently, a fresh interest in the topic has been rekindled by the work of 
Morris and Thorne |@, leading to a flurry of activity branching off into diverse directions. A 
brief resume of current work devoted to the physics of Minkowski-signature wormholes includes 
topics addressing fundamental features of traversable wormholes iS, explicit modelling of 
wormhole metrics and the corresponding classical P] and quantum mechanical stability [0 
analyses, wormholes as time machines and the problem of causality violation |^, wormholes in 
higher-derivative gravity [^], wormholes from the gravitationally squeezed vacuum []^, possible 
cosmological consequences of early universe wormholes H.Q. and wormholes as gravitational 


lenses |]R]. A thourough and up-to-date survey of the present status of Lorentzian wormholes 
may be found in the excellent monograph by Visser . 

There are plausible physical arguments suggesting that Lorentzian wormholes should exist 
at least at scales of order the Planck length. Most of what is known about them is based on 
detailed analyses of models, and within the literature devoted to the subject, the existence of 
wormholes is taken as a working hypothesis. Metrics describing wormholes with desirable traits 
are written down by hat and the properties of the corresponding hypothetical stress-energy 
tensors needed to support the wormhole spacetime are then worked out and analyzed. In an 


example of an analysis of this sort. Ford and Roman |]^ have derived approximate constraints 
on the magnitude and duration of the negative energy densities which must be observed by a 
timelike geodesic observer in static spherically symmetric wormhole spacetimes. More recently, 
Taylor, Hiscock and Anderson have argued that stress tensors for massive minimally and/or 
conformally coupled scalars fail to meet the requirements for maintaining hve particular types 
of static spherically symmetric wormholes, but have not solved the backreaction problem [^. 
In particular, no one up to now has succeeded in writing down a bona-hde wormhole solution of 
either the classical or semiclassical Einstein held equations. The reason for this state of affairs is 
easy to understand. In the hrst case, it is well known that any stress-energy that might give rise 
to a wormhole must violate one or more of the cherished energy conditions of classical general 
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relativity ii- Hence wormholes cannot arise as solutions of classical relativity and matter. 
If they exist, they must belong to the realm of semiclassical or perhaps quantum gravity. In 
the realm of semiclassical gravity, one sets the Einstein tensor equal to the expectation value 
of the stress-energy tensor operator of the quantized helds present: 




( 1 ) 


A primary technical difficulty in semiclassical gravity is that (T^y) depends strongly on the 
metric and is generally difficult to calculate. Until recently, all calculations of have been 
performed on hxed classical backgrounds. The hxed background in turn, as its name implies, 
must be a solution of the classical Einstein equation. As there are no classical wormhole 
backgrounds, no corresponding semiclassical back-reaction problem can be set up meaningfully. 


In this Letter we present and summarize the results of the hrst self-consistent wormhole 
solutions of semiclassical Einstein gravity. Prior to this, a self-consistent wormhole solution 
had been obtained using a phenomenological stress tensor not derived from quantum held 
theory [j^. The results of the present calculation may be taken as numerical evidence for the 
existence of Lorentzian wormholes. For the source term in (1) we employ the stress-energy 
tensor of Anderson, Hiscock and Samuel, which is calculated for a quantized scalar held in 


an arbitrary static and spherically symmetric spacetime |^. This means that in the held 
equation (1), both the Einstein tensor as well as the stress tensor individually depend on two 
independent functions of the radial coordinate. When supplemented with the appropriate set 
of boundary conditions, the solutions of the resultant coupled nonlinear diherential equations 
are therefore self-consistent because both the spacetime metric and the distribution of stress- 
energy are determined simultaneously and coherently. This should be contrasted clearly with, 
and distinguished from, the approach taken in perturbative back-reaction problems in which 
the (background) spacetime is hxed once and for all, and the stress tensor is supplied as an 
explicitly known function of that hxed background. 


In the following, we set up the semiclassical held equations valid for any static and spherically 
symmetric spacetime containing quantized scalar matter and discuss the nature of the boundary 
conditions needed for solving this system of fourth-order equations. We consider the case of a 
conformally coupled scalar held. Results of the numerical calculations are presented graphically 


3 










and we appeal to an approximate but analytic treatment to reveal the important asymptotic 
behavior of the solution. Further details of this and related calculations will appear in a separate 
paper. 


The metric for a general static and spherically symmetric spacetime can be cast into the 
form 

= —f{l)dt‘^ + df + r^(/) (dO"^ + sin^ 9 dcj)'^^ , (2) 

where f{l),r{l) are two independent functions of the proper distance 1. This form of the metric 
is suitable for dealing with wormholes, or for that matter, any static and spherically symmetric 
spacetime which might contain a throat, i.e., r(0) > 0. Following the established nomenclature, 
/ is denoted the redshift function and r is known as the shape function |]T3. In this metric. 


the semiclassical field equations take the form HlQ 


G;:(/(0,r(0) = 87r(r;(/(/),r(0;6), (3) 

where ^ is the (non-minimal) scalar coupling to the metric. Note the dependence of both sides 
of the equation (3) on the two unknown functions. In the metric (2), the components of the 
Einstein tensor are given by Gj = ^ ^ ^ and + ^ + 

the prime denotes the derivative with respect to 1. An accurate analytic approximation to the 


exact numerically calculated scalar field stress energy tensor was developed in Ref. and is 
expressed there as 

= (t;)„ + (e - i)(r;), + ((- (4) 

where the individual factors (T’^)o,i ,2 are written in terms of two functions of the radial co¬ 
ordinate. The last factor (T^)zo 3 was left in terms of combinations of curvature tensors and 
covariant derivatives |^. To be of practical use in the present calculation, we must work out 
these curvature terms and transform all the factors in (4) in terms of our two functions /(/), r(/) 
of the proper distance 1. Carrying out these straightforward but rather lengthy steps we find 
that the semiclassical Einstein equations for a conformally coupled scalar (^ = ^) are as follows 

('R'2 ^ 1 )■ 

^ 57607r^- 


the tt component: 


2r" 

r 
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the 69 component: 
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The full expressions for the general ^-dependent components of the Einstein equations are 
too lengthy to be shown here. 

In order to solve the above held equations, we must supply an appropriate set of bound¬ 
ary conditions. Note that the tt and 6'0-equations are fourth order differential equations in 


5 



the two functions / and r while the //-equation is third order. This latter equation is actu¬ 
ally a constraint which plays the role of restricting the solutions of the coupled fourth order 
differential equations |^. Since we seek wormhole solutions of (1), we shall specify the bound¬ 
ary data at / = 0, the origin of proper distance as measured from the throat. The complete 
set of boundary conditions therefore requires specifying eight pieces of information, namely 
/(O),/'(O),/"(O),/'"(O) as well as r(0), r'(0), r"(0) and r'"(0). In particular, r(0) is the (worm- 
hole) throat radius, and we can use the //-equation to derive an exact relation between r(0) and 
just three initial conditions. In order to do so, let us recall the boundary conditions appropriate 
to a wormhole (or more generally, for any spacetime with a throat). If a solution of (1) is to 
possess a throat, then we must require that r(0) > 0, r'(0) = 0 and r"(0) > 0. In other words, 
there must exist a sphere of minimum radius located at the origin of proper distance. So r(/) 
is simply a positive increasing function of / in the neighborhood of the origin. For two-way 
passage through the wormhole, it is judicious to avoid solutions with horizons at the throat. 
This is taken care of by requiring a (locally) nonvanishing redshift function: /(O) > 0 (the other 
possibility, /(O) < 0, leads to a Euclidean metric). The chain rule also provides an additional 
constraint: namely r'(0) = 0 — f'{0) = 0 0. We may ask that the redshift be a locally 
increasing function: f'\0) > 0. These constitute the minimum requirements. Beyond this, one 
could impose a symmetry on the solutions of the form r(/) = r(—/) and /(/) = /(—/) which 
automatically eliminates all the odd derivatives: = 0. In particular, one 

can set r'"(0) = f"'{0) = 0. Taken together, these give the boundary conditions appropriate 
to a wormhole. However, one can not choose freely all the eight boundary conditions indepen¬ 
dently. This is easily seen by writing out the //-equation (6) at the point / = 0, which yields an 
algebraic quartic equation for the throat radius r(0): 


-4(^)'[i+i.i(/(o))Hor 


+ 32 


|i + bn/(o))Mo)= 


( 8 ) 


+ [R-^ - 16r"(0)" ln(/(0))]r(0)" + 16 hi(/(0)) = 0, 


where we have used only that r'(0) = f'{0) = 0. Indeed, it is important to note that this 
constraining relation is independent of the values assigned to the third derivatives f''\0) and 
r'"(0). With some effort, this quartic equation can be solved for r(0), though the resulting 
expressions are not particularly transparent. However, it suffices for our purposes to consider 


6 





(8) in certain simplifying but natural cases in order to get a feeling for the allowed range in 
r(0). In this way, we have found that large throat radii can result even for values of /(O), f'\0) 
and r"(0) of order unity. As case in point, taking /(O) = f"{0) = 1, r"(0) = 0 and = | yields 
r(0) ~ 67/p where Ip is the Planck length. Larger (and smaller) throat radii are also possible 
depending on the values chosen for the other boundary data and the scalar coupling constant 

e 

We now turn to our specihc calculations of r(/) and /(/) subject to the above class of 
boundary conditions. Employing the Rugge-Kutta method, we have obtained a good numerical 
solution of these equations using the following boundary conditions for /(/) and r(/) at the 
throat I = 0: 

-l<ln/(0)<0; /'(0) = 0; /"(O) = 0; /'" (0) = 0; (9) 

r(0) = ^-16772 ln/(0); r'(0) = 0; r"(0) = 0; r'"(0) = 0 

The graphs which represent a particular numerical solution for ln/(0) = —| for various length 
scales are displayed below. The redshift function /(/) is plotted versus proper distance in 
Figures la and lb; both axes have been scaled by K~^. In Fig. lb, we see that /(/) is 
a positive increasing function. Actually, there is a hne-structure in the form of Planck-scale 
“ripples” or spatial oscillations superimposed on this function; these are not numerical artifacts. 
These oscillations are clearly revealed in the “blow-up” graph in Fig. la. The wormhole’s 
shape function r(/) is depicted in Figs. 2a and 2b. In Fig. 2a one sees the throat, of radius 
r(0) ~ 0.02 Ip, in the neighborhood of the origin extending out to about lOA'"^ at which point 
the wormhole “flares” out, marking the onset of the superimposed small-scale oscillations. The 
graph in Fig. 2b shows the gross features of the shape function. As remarked below, the 
oscillations in both /(/) and r{l) are composed of two different modes. Of course, since both 
/ and r are even functions, their graphs can be reflected through the origin: I —> —1. The 
wormhole’s embedding diagram is easily inferred from the Figures 2a and 2b. 

It is important to know the asymptotic behavior of /(/) and r(/). To get at this information, 
we carried out an asymptotic analysis of the Einstein equations (5-7) taking into account and 
guided by the results of the numerical investigation. From the numerical calculations, we see 
that for sufficiently large /, both the redshift and shape functions can be represented as a sum 
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of two distinct components /(/) = F{1) + 0(/); r(/) = R{1) + p(/), where F{1) and R{1) are 

strictly monotone increasing and ^(/) and p{l) are bounded oscillating functions. 

The relative magnitudes of these components and their derivatives may be estimated straight¬ 
forwardly and then used to expand consistently the coupled Einstein equations. We hnd 
that the oscillating modulation is composed of two modes with frequencies uf = and 
'^2 = 16 X 2 ( 4 + 3 InF) ’ I’Gspectively, while in the limit of large I, R{1) ~ I and F{1) ^ (a In / — 6)^ 
for a = 5.3 and b = 25.5. 

From these combined numeric and analytic calculations we see that for the chosen set of 
boundary conditions, R{1) I, thus our self-consistent wormhole connects two spatial regions 
which are asymptotically flat, modulo the Planck-scale wiggles. The redshift function however, 
does not approach a constant value as / —> cx), so the metric as a whole is not asymptotically 
flat. We have found additional self-consistent solutions of (1) by taking different values for 
the scalar coupling and the boundary data. In this way, we have found local solutions which 
correspond to large throat (r(0) ~ 200 — 300 Ip) wormholes with horizons located far from 
the throat and wormholes connecting two bounded spatial regions. A full account of these 
calculations will appear in a separate publication. 

The work of S.V.S and A.P. was supported in part by the Russian Foundation of Funda¬ 
mental Research grant No. 96-02-17366a. 
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Figure la: The redshift function on small scales 
Figure lb; The redshift function on large scales 


11 



Figure 2a: The wormhole shape function on small scales 
Figure 2b; The wormhole shape function on large scales 
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